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ABSTRACT 

We  develop  a  method  for  the  analysis  of  perturbations  of  integrable 
planar  systems  of  differential  equations.  Concentrating  on  the  case  in 
which  the  unperturbed  system  is  Hamiltonian  and  the  perturbation  introduces 
dissipation  and  time-periodic  forcing,  we  show  how  the  global  solution 
curves  (level  sets)  of  the  unperturbed  system  can  be  used  in  regular  por- 
turbation  calculations  to  locate  subharmonic  orbits  and  homoclinic  orbits 
and  to  characterize  the  bifurcations  in  which  they  are  created  as  external 
parameters  are  varied.  We  apply  our  results  to  Duffing's  equation  and 
point  out  applications  to  the  chaotic  motions  of  buckled  elastic  beams 
undergoing  periodic  excitation. 
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1.  INTRODUCTION:  A  PHYSICAL  EXAMPLE 


In  this  article  we  develop  a  method  for  the  analysis  of  systems  of  the 

form 

x  =  f(x)  +  eg(x;t),  (1.1) 

r  'i 

fu'i  2  fi(u,v)  9i  (u,v;t) 

where  x  -  [VJ  e  IR  ,  and  f  =  .  9  “  [g2(u,vit)j  are  suf«'iently 

smooth  (bounded)  functions  and  g  is  T-periodic  in  t.  In  examples  of  in¬ 
terest  g  will  generally  also  contain  parameters 

g^(Xjt)  =  g(x;t;u)  u  e  IRm  .  (1 .2) 

The  unperturbed  system, 

x  =  f ( x ) »  (1.3) 

is  assumed  to  be  integrable.  In  particular,  we  shall  concentrate  on  the 

2 

Hamiltonian  case,  in  which  a  real  valued  function  H  :  IR  e  IR  exists  and 

f,(u,v)=f,  f2(U.*)  ■  -  ft  •  0.4) 

The  results  we  describe  here  can  be  generalized  to  multidimensional 
systems  and  even  to  infinite  dimensional  Hamiltonian  systems  arising  on 
evolution  equations  in  studies  of  partial  differential  equations  (cf.  Holmes 
and  Marsden  [1981a]),  and  to  the  study  of  multidegree  of  freedom  Hamiltonian 
systems  with  perturbations  which  do  not  depend  explicitly  on  time  (cf. 
Churchhill  [1980],  Holmes  and  Marsden  [1981b, c]). 

Our  basic  assumption,  made  precise  in  Section  3,  is  that  the  unper¬ 
turbed  system  (1.3)  possesses  a  nondegenerate  (hyperbolic)  saddle  point 
pQ  wit.i  a  homoclinic  orbit  T°  =  (g°(t)|t  6  IR)  such  that 


1.2 


- 1 


« 


♦ 


lim  q°(t)  =  lim  q°(t)  =  p  .  (1.5) 

t-+  +«  t**-  ”®  ® 

We  further  assume  that  the  interior  of  r°  is  filled  with  a  one-parameter 
family  of  periodic  orbits  ra  =  {qa(t)|a  €  (-1,0),  t  €  [O.T^]}  where  Tq 

ct 

is  the  period  of  q  (t)> 

lim  qa(t)  =  q°(t) ,  (1.6) 

a->-  0 

a  -1 

and  lim  q  (t)  =  q  is  an  elliptic  fixed  point  or  center,  cf.  Figure  1.1. 
a-*-- 1 

The  homoclinic  orbit  can  be  replaced  by  a  homoclinic  cycle  connecting  several 


qa(t) 


Figure  1.1.  The  unperturbed  phase  portrait  of  (1.3). 

saddle  points.  Such  oroits  and  cycles  filled  with  periodic  motions  occur 
naturally  in  the  phase  planes  of  single  degree  of  freedom  nonlinear  Hamil¬ 
tonian  systems;  the  pendulum: 

G  *  sin  G  °  0  ,  (1.7) 

and  Ouffing’s  equation  with  negatv.a  linear  stiffness: 

y  -  y  ♦  y3  ■  0,  (1.8) 

provide  examples.  We  select  Duffing' s  equation  as  our  major  example  in 
subsequent  sections  because  of  its  importance  in  mechanics,  which  we  now 


1.3 


briefly  indicate.  For  more  details  see  Tseng  and  Dugundjii  [1971],  Moon 
and  Holmes  [1979],  Holmes  [1979],  Moon  [1980])  and  Holmes  and  Marsden 
[1981a]. 

Consider  a  slender  elastic  rod,  clamped  at  one  end  in  a  rigid  frame¬ 
work  and  constrained  to  move  in  a  plane  (Figure  1.2).  Two  magnets  attached 
to  the  frame  as  indicated  cause  the  beam  to  buckle  either  to  the  left  or 


figure  1.2.  The  maqnetoelastic  beam 

(from  Moon  and  Holmes  [1979]) 


right,  the  central  position  being  unstable.  The  whole  framework  is  now 
moved  sinusoidally,  so  that  the  beam  vibrates  under  'inertial'  excitation, 
for  sr.«al!  excitation  amplitudes,  periodic  motions  are  observed  about  either 
stable  equilibrium,  but  as  the  assplitude  is  increased  an  apparently  sudden 
iratsition  to  a  ‘chaotic  snap  through’  motion  is  observed,  in  this  state, 
the  bear,  oscillates  irregularly  about  first  one  and  then  the  other  equilib¬ 
rium.  An  exaroic  of  such  a  motion, measured  from  a  strain  guage  at  the 


beam's  root,  is  shown  in  Figure  1.3. 


Figure  1.3.  Nonperiodic  motion  of  (a)  the  beam  (from  Moon  and 
Holmes  [1979])  and  (b)  the  forced  Duffing  equation 
(from  Holmes  [1979]) . 

As  shown  by  Moon  and  Holmes  [1979],  the  nonlinear  partial  differential 
equations  of  an  inextensible  elastic  rod  in  a  nonuniform  magnetic  field 
can  be  truncated  in  a  single  mode  Galerkin  approximation  to  yield  Duffing's 
equation: 

y  _  y  +  y3  s  0  ,  (1.9) 

wnere  y  *  y(t)  is  the  (nondimensional )  amplitude  of  the  first  mode  of 
vibration  and  dissipation  is  ignored.  On  the  introduction  of  weak  dissipa¬ 
tion  (aerodynamic  damping)  and  periodic  forcing  one  obtains 

y  -  y  f  y3  =  £(-f  COS  wt  -  )  .  (1.10) 

which  is  tne  equation  to  be  studied  in  Sections  4-6  of  this  article- 

in  the  following  section  (2)  we  introduce  some  b*sic  concepts  from  the 
qualitative  theory  of  dynamical  systems  and  then  in  Section  3  state  and 
sketch  the  proofs  of  the  results  on  perturbed  intecrable  systems.  In  sec¬ 
tion  4  »,e  app  1  y  these  results  to  Duffing's  equation  and  then,  in  Section  5, 


1.5 


review  more  of  the  relevant  abstract  theory.  Finally,  in  Section  6,  we 
bring  the  abstract  theory  back  to  bear  on  Duffing's  equation  in  an  attempt 
to  understand  irregular  motions  such  as  that  of  Figure  1.3.  In  this  article, 
although  we  assume  some  familiarity  with  the  qualitative  theory  of  dynami¬ 
cal  systems  and  bifurcation  theory,  we  do  attempt  to  introduce  or  review  most 
of  the  important  aspects  necessary  for  our  study.  Background  material  on 
differential  dynamics  may  be  found  in  Chillingworth  [1976],  and  on  bifurca¬ 
tion  theory  and  nonlinear  analysis  in  Iooss  and  Joseph  [1980]  and  Chow  and 
Hale  [1981].  The  latter  reference,  and  Chow,  Hale  and  Mallet-Paret  [1980] 
contain  many  of  the  results  of  our  Section  3,  but  presented  from  a  different 
viewpoint.  The  book  of  Andronov,  Vittand  Khaiken  [1966]  continues  to  pro¬ 
vide  the  best  background  in  nonlinear  oscillations. 


2.1 


2.  POINCARE  MAPS  AND  INVARIANT  MANIFOLDS 


We  consider  the  periodically  forced  system 


x  =  f(x)  +  eg(x,t);  x  S  IR2 ,  glx.t)  =  g(x,t  +  T)  ,  (2.1) 

which  can  be  rewritten  as  an  autonomous  system  defined  on  the  'toroidal'  phase 

2  1 
space  IR  x  S  : 


x  =  f(x)  +  eg(xt0) ' 
9  =  1 


(x,e)  e  ir2 


(2.2) 


where  =  IR/T,  the  circle  of  length  T.  This  naturally  copes  with  the 
T-periodicity  of  the  vector  field.  We  next  define  a  cross  section  for  the 
flow: 

t 

i  =  {(x,e)  €  IR2  x  s1 le  =  t0  e  [o,T)},  (2.3) 


and  a  first  return  or-  Poincarl  map 


e 


obtained  by  following  orbits  (x(t+  t  ),  9(t  +  t  )) 
t  oo 

(x(t  ),  tQ)  6  E  0  to  their  next  intersection  with 


based  at  points 
t 

Z  °.  We  have 


(2.4) 


t 

P£°(x(0))  =  ir* (x(T  +  tQ),  T  +  tQ)  (2.5) 


where  (x(t  +  tQ),  t  +  tQ)  is  the  solution  of  (2.2)  based  at  (x(tQ),  t  ) 

and  tt  denotes  projection  onto  the  first  factor.  When  e  =  0,  the  unper- 

t  t 

turbed  Poincare  map  PQ  =  PQ°  is  identical  on  every  section  Z  0  in  view 
of  the  invariance  of  x  =  f(x)  under  time  translations.  PQ  is  given  simply 
by  the  time  T  flow  map  of  the  unperturbed  equation 


2.2 


Pp(x(o))  «  x°(T)  ,  (2.6) 

where  x°(t)  is  the  unperturbed  solution  of  2.1  based  at  x(0).  The  per- 
t 

turbed  maps  Pe°  differ  on  different  sections  but  any  two  are  diffeomorphic 

t  +nT  t 

and  P  0  =  P  o. 

e  e 

The  Poircci'd  map  captures  the  dynamics  of  (2.1)  in  the  sense  that  T- 

t 

periodic  orbits  of  (2.1)  correspond  to  fixed  points  of  P£°,  mT-periodic 

subharmonics  correspond  to  periodic  cycles  of  period  m  and  stability  types 

correspond.  The  study  of  fixed  and  periodic  points  of  the  two  dimensional 
+■ 

map  P£°  is  generally  uasier  than  that  of  the  corresponding  periodic  motions 
in  the  three  dimensional  flow  of  (2.2). 

Suppose  that  the  map  P  (the  sub-and  superscripts  will  be  dropped 
where  they  are  not  explicitly  needed)  has  a  hyperbolic  fixed  point  p.  This 
implies  that  the  linearized  map  OP(p)  has  eigenvalues  X i  with  modulus 
|X^  |  M *  i  =  1 ♦  2:  a  sink,  a  saddle  or  a  source.  The  stable  manifold 
theorem  (cf.  Chill ingworth  [1976],  Hartman  [1973],  Hi rsh-Pugh-Shub  [1977]) 
and  Appendix  A)  asserts  that,  in  a  neighborhood  U(p)  there  are  smooth  local 
stable  manifolds  W^QC(p),  wlj1  (p)  tangent  to  the  eigenspaces  Es,  Eu  of 
the  linearized  problem 

i  =  DP(pK  ,  (2.7) 

and  of  dimension  s,  u  (s  +  u  *  2  here)  where  s(  resp.  u)  is  the  number  of 
eigenvalues  with  moduli  less  than  1  (resp.  greater  than  i).  These  manifolds 
are  simply  defined  as  the  sets  of  points  asymptotic  to  p  under  forward 
(resp.  backward)  iterations  of  P 

W^QC(p)  =  (x  e  U(p) |Pn(x)  e  u(p),  Vn,  and  Pn(x)  ■*  p  as  n  -  «}, 
wJ0C(p)  -  {X  €U(p)jp-n(x)  €U(p),vn,  and  P"n(x)  -  p  as  n  -*■«}, 


and  global  manifolds  are  obtained  by  iteration  of  the  local  manifolds 
backwards  (resp.  forwards): 


NS(P)  -  u  P'V  (p))  . 
n>0  loc 


w“(p)  •  u  P>“  (p)) 

n>0  ,oc 


Of  course  if  both  eigenvalues  lie  inside  (resp.  outside)  the  unit  circle 
then  Wu(p)  (resp.  Ws(p))  is  empty. 

We  note  that  while  locally  they  are  smooth  submanifolds  of  IR^, 

S  ii 

W  (p)  and  W  (p)  are  not  necessarily  submanifolds  globally,  since  they 
nay  wird  back  and  forth  and  accumulate  on  themselves,  as  they  do  in  the 
specific  examples  to  follow  (Figure  2.1). 


It  is  important  to  realize  that  the  orbit  of  a  map  (Pn(x))  is  a 
sequence  of  points ,  rather  than  a  curve  *s  in  the  case  of  orbits  of  flows 
arising  from  vector  fields.  Each  of  the  curves  Ws(p',  Wu(p)  in  Figure 
2.1  contains  a  one  parameter  family  of  such  orbits.  Uniqueness  of  solutions 


implies  that  the  stable  (or  unstable)  manifolds  Ws(p^),  W5^)  of  two 
distinct  points  cannot  intersect  and  nor  can  such  a  manifold  intersect 
itself. 

A  homoclinic  point  q  to  a  hyperbolic  saddle  point  p  is  a  point 
whose  orbit  approaches  p  under  both  forward  and  backward  iterates  of  P: 

lim  Pn(q)  s  lim  P~n{q)  3  p, 

n* »  tv*  * 

thus  q  €  wu(p)  nws(p).  When  two  plane  curves  intersect  at  a  point  they 
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generally  do  so  transversely,  so  that  their  tangent  vectors  span  IR  . 

Such  a  point  q  €  Wu(p)  ^Ws(p)  is  called  a  transversal  homoclinic  point 
and  will  play  an  important  role  in  our  analysis.  The  existence  of  one  such 
point  immediately  implies  the  existence  of  infinitely  many,  lying  on  the 
orbit  {Pn(q)L  Moreover,  the  X-lemma  (Palis  [1969],  Newhouse  [1980]  and 
Appendix  A)implies  that  W  (p)  accumulates  upon  itself,  and  W  (p)  accu¬ 
mulates  upon  itself,  leading  to  the  complicated  structure  of  Figure  2.2. 


Figure  2.2.  Transverse  homoclinic  orbits.  See  Section  5.1 
for  discussion  of  the  rdle  of  the  rectangle  S. 


In  Section  5  we  discuss  some  implications  of  such  homoclinic  orbits, 
but  first  we  develop  a  method  for  their  detection  in  specific  examples. 
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3.  PERTURBATIONS  OF  INTEGRABLE  SYSTEMS:  MELNIKOV'S  METHOD 


The  basic  ideas  to  be  developed  in  this  section  are  due  to  Melnikov 
[1963],  More  recently  Chow,  Hale  and  Mallet-Paret  [1980]  have  obtained 
similar  results  using  alternative  methods,  and  Holmes  and  Marsden  [1981a, b,c] 
have  applied  the  method  to  certain  infinite  dimensional  flows  arising  from 
partial  differential  equations  and  to  multidegree  of  freedom  autonomous 
Hamiltonian  systems.  The  basic  idea  is  to  make  use  of  the  globally  compu¬ 
table  solutions  of  the  unperturbed  inteqrable  system  in  the  computation  of 
perturbed  solutions.  To  do  this  we  must  first  ensure  that  the  perturbation 
calculations  are  uniformly  valid  on  arbitrarily  long  or  semi-infinite  time 
intervals. 

First  we  make  our  assumptions  precise.  We  consider  systems  of  the  form 


e  IR 


(3.1) 


where  f  = 


V*)] 

f2(x) 


g  = 


x  =  f (x)  +  eg(x,t);  x  = 

9i(x»t)l 

g  (x  t)J  are  su-,:^ic1'ent1Y  smooth  (C  ,  r  _>  2)  and 


bounded  on  bounded  sets  and  g  is  T-periodic  in  t.  For  simplicity  we  assume 


-  3H 


3H 


that  the  unperturbed  system  is  Hamiltonian  with  f-j  =  ^2  =  “  fu'  ”^e 

non-Hamiltonian  case  is  considered  by  Melnikov  [1963]  and  Holmes  [1980]. 
Specific  assumptions  on  the  unperturbed  flow  are: 


A1 .  For  e  =  0  (3.1)  possesses  a  homoclinic  orbit  q°(t)  to  a  hyperbolic 
saddle  point  p  . 

A2.  Let  r°  =  (q°(t)|t  €  IR }  u  {pQ}.  The  interior  of  r°  is  filled  with 

a  continuous  family  of  periodic  orbits  qa(t),  a  s  (-1,0).  Letting 

d(x,r)  =  inf  { x-q |  we  have  lim  sup  d(qa(t),  r)  =  0. 
q  Sf  a-*-  0  t£  IR 
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A3.  Let  h^  =  H(qa( t) )  and  be  the  period  of  qa(t).  Then  Ta  is  a 

differentiable  function  of  h  and  dT  /dh  >  0  inside  r°. 

a  a  a 

mono  tcflitfdty, 

We  note  that  A2  and  A3  imply  that  -»•  ®^as  a  0.  Many  of 
the  results  to  follow  can  be  proved  under  less  restrictive  assumptions.  In 
what  follows  we  indicate  the  proofs  but  omit  some  details.  See  Greenspan 
[1981]  for  more  information. 

3.1  Bifurcations  to  Homoclinic  Orbits 

Lemma  3.1 .  Under  the  above  assumptions,  for  e  sufficiently  small  (3.1) 

has  a  unique  hyperbolic  periodic  orbit  y°(t)  =  p  +  0(e).  Correspond! nqly, 

t  to 

the  Poincarg  map  P£°  has  a  unique  hyperbolic  saddle  point  P£°  =  PQ  +  0(e). 

Proof.  This  is  a  straightforward  application  of  the  implicit  function 
theorem,  our  assumptions  implying  that  DPq(p0)  does  not  contain  1  in  its 
spectrum  and  hence  that  Id  -  DPq(pq)  is  invertible  and  there  is  a  smooth 
curve  of  fixed  points  in  (x,e)  space  passing  through  (pQ,0).  ■ 

Lemma  3.2.  The  local  stable  and  unstable  manifolds  W?  (v  ) ,  wV  (y  ) 
-  -  loc^e'’  loc  e 

of  the  perturbed  periodic  orbit  are  Cr  close  to  those  of  the  unperturbed 
periodic  orbit  pQ  x  s\  Moreover,  orbits  q£(t,tQ),  q£(t,tQ)  lying  in 
VIs  (yJ,  Wu(Ye)  and  based  on  E  0  can  be  expressed  as  follows,  with  uniform 
validity  in  tne  indicated  time  intervals. 


q*(t,t0)  »  q°(t  -  t0)  *  cq*(t,t0)  ♦  G(e2),  t  €  [tQ»»); 
q“(t,t0)  »  q°(t  -  t0)  +  cqJ(t,tQ)  +  0(e2),  t  6  (-®,tQ]. 


Proof.  The  existence  of  the  perturbed  manifolds  follows  from  the  invariant 


manifold  theory  of  Hirsch,  Pugh  and  Shub  [1977],  cf.  Hartman  [1973].  We 
fix  a  v-neighborhood  (0  <  e  «  v  «1)  Uv  of  pQ  inside  which  the  local 
perturbed  manifolds  and  their  tangent  spaces  are  e-close  to  those  of  the 
unperturbed  flow  (or  Poincard  map).  A  standard  Gronwall  estimate  shows  that 
perturbed  orbits  starting  within  0(e)  of  q°(o)  remain  within  0(e)  of 
q°(t  -  t  )  for  finite  times  and  hence  that  one  can  follow  any  such  orbit 
from  an  arbitrary  point  q°(0)  on  r°  outside  Uy  to  the  boundary  of  Uy, 
at,  say  t  *  t-j .  Once  in  Uy,  if  the  perturbed  orbit  q£  is  selected  to 
lie  in  Wb(Y£),  then  its  behavior  is  governed  by  the  exponential  contrac¬ 
tion  associated  with  the  linearized  system.  Straightforward  estimates  then 
show  that  |qS(t,tQ)  -  q°(t  -  t  )|  =  0(e)  for  t  €  (t|,«).  Reversing  time, 
one  obtains  a  similar  result  for  q^( t , t  ) .  ■ 

Sanders  [1980] vas  the  first  to  work  out  the  asymptotics  in  detail. 

Tnis  lemma  implies  that  solutions  lying  in  the  stable  manifold  are  uniformly 
approximated  by  the  solution  q^  of  the  first  variational  equation: 

q*(t,tQ)  »  Df(q°(t  -  tQ))q*(t,t0)  +  g(q°(t  -  tQ) ,t)  .  (3.3) 

A  similar  expression  holds  for  q^1  ( t,tQ) .  Note  that  the  initial  time, 
t  appears  explicitly  since  solutions  of  the  perturbed  systems  are  not 
invariant  under  arbitrary  translations  in  time. 

t  t 

We  next  define  the  separation  of  the  manifolds  Wu(p  °),  Ws(p  ) 

^  £  £ 

on  the  section  z  0  at  the  point  q°(0)  as 

d(tQ)  °  q“(t0)  -  q|(tQ)  ,  (3.4) 

where  q“(tQ)  «  ^VV*qe*V  ^  q*(t0,tQ)  are  the  unique  points 


t  t  t 

on  Wu(p°),  Ws(p,°)  'closest'  to  P  0  and  lying  on  the  normal 
fVeo))  =  (f2(q°(o)),  (q° (o ) ) )T  to  r  at  q°(o).  The  Cr  closeness 

of  the  manifolds  to  r,  and  Lemma  2  then  imply  that 


d(tQ)  =  e 


f(q°> :  ♦  0  (e2) 


(3.5) 


Here  the  wedge  product  is  defined  a  »  b  s  •  agb-,  and  f  *  (q?  -  q^) 
is  the  projection  of  q^j1  -  q*  onto  f\  cf.  Figure  3.1. 


«“<p*°) 


>^?po 

".C\x 


w  - 

s  vW 


\  q°(0)  *  eqj(t  ) 

\\  /  ° 

A/7 

\  Aqu(tn) 

\c  Ve  o'  I 

— 1-4- - (q°(0)) 


q°(0)  +eq^(t0) 


wxVx: 

e  e  X 


^(q°(0)) 


Figure  3.1  The  perturbed  manifolds  and  the  distance  function. 


Finally  we  define  the  Melnikov  function 


M(t  )  *  P  f(q°(t  -  tQ))  -  g(q°(t  -  tQ),t)  dt  .  (3.6) 

'  *03 


Theorem  3.3.  If  M(t  )  has  simple  zeros  and  maxima  and  minima  of  0(1),  then 

_  o  t  i 

for  z  >  o  sufficiently  small  Wu(p°)  and  Ws(p^°)  intersect  transversely. 

£  t 

l_f  M(t^)  remains  bounded  away  from  zero  then  Wu(p£o)  nws(p£°)  3  $. 
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Proof.  Consider  the  time  dependent  distance  function 


Mt,t0)  *  f(q°(t  -  tQ))  -  (qf(t,t0)  -  qf(t,t0)  (3.7) 


AU(t,t0)-AS(t,t0) 


and  note  that  d(t  )  3  eA(t  ,t  )/|f(q°(*j))  j  +  G(c~),  from  (3.5).  We  com- 


AS(t,t0)  =  Of (q°(t  -  t0))qQ(t  -  tQ)  -  q*(t,tQ) 


+  f(q°(t  -  t  ))  -  q^(t.t  ) 


Using  (3.3)  and  the  fact  that  qQ  =  f(qQ)  this  yields 


As  •  Df(q°)f{q°)  -  q*  +  f(q°)  -  (Of(q")qJ  +  g(q\t)) 
=  trace  Of (q° )AS  +  f(q°)  *  g(q°,t) 


Os.S  .  /  0 


(3.8) 


But,  since  f  is  Hamiltonian,  trace  Of  =  0.  Integrating  (3.8)  from  t 

V 

to  we  have 


ASKtQ)  -  As(t0,tQ)  =  f(q°(t  -  tQ))  -  g(q°(t  -  tQ),t)  dt.  (3.9) 

to 


However  AS(®,t  )  =  lim  f(q°(t  -  t  ))  *  qf(«,t  )  and  lim  q°(t-t  )  9  p  , 

0  £  aa  ®  10  %  +  a>  0  ^ 

so  that  lim  f(q°(t  -  t.))  3  0  while  qf (°°. )  is  bounded,  from  Lemma 

t-*  » 

3.2-  thus  AS(t“, tQ)  °  0  and  (3.9)  gives  us  a  formula  for  AS(t0,tQ).  A 
similar  calculation  gives 
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=  i  f(q0(t  ■  v5  a  g(q0(t  ■  t0)*t)  dt*  (3-'°) 
'  *00 


and  addition  of  (3.9)  and  (3.10)  and  use  of  (3.5)  yields 


rf(tQ)  «  cM(to)/|f(q°(0))|  +  O(e-)  . 


(3.11) 


Since  |f(q°(0))|  =  0(1),  M(t  )  provides  a  good  measure  of  the 

0  t 

separation  of  the  manifolds  at  q°(0)  on  Z  In  particular,  if 

M(tQ)  oscillates  about  zero  with  maxima  and  minima  independent  of  e, 

then  q“(t  )  and  q*(t  )  must  change  their  orieotiation  with  respect 

to  (q°( 0) )  as  t°  varies  and  hence  Wu(Pe°)  and  WS(P£°)  intersect. 

~~  t  o] ,  then  it  follows  that  the  intersections 
<dto  > 

are  transversal.  Conversely,  if  no  zeros  exist,  then  q^  and  q*  retain 


If  the  zeros  are  simple 


the  same  orientation  and  hence  the  manifolds  do  not  intersect. 


Remarks  1.  We  note  that  M(t  )  is  T-pericdic  in  t. ,  as  it  should  be. 
-  o  o 

t  t  +T 

since  the  maps  P_°  and  P£0  are  identical  and  thus  d(tQ)  *  d(tQ  +  T). 
In  computing  M(t  )  we  are  effectively  stana.ng  at  i  fixed  point  q°(0) 
on  a  moving  cross  section  Z  0  and  watching  the  perturbed  manifolds  oscil¬ 
late  as  t  varies.  In  his  analysis,  Greenspan  [1981]  fixes  the  section 
and  moves  the  base  point  q°(0). 

2.  If  the  perturbation  g  is  derived  from  a  (time  dependent)  Hami1- 
tonian  function  G,  g1  =  — g^  =  -  then  we  have 

M(t0)  »  -f  (H(q°(t  -  tQ)).  G(q°(t  -  t  ),t))  dt  ,  (3.12) 

where  ( • ,  •}  denotes  the' Poisson  bracket.  This  is  a  natural  formula  if 
one  recalls  that  the  first  variation  of  the  unperturbed  1  ton i an  energy 


3.7 


H  will  be  obtained  by  integrating  its  evolution  equation 


H  =  {H,G} 


(3.13) 


along  the  unperturbed  orbit  q°(t  -  tQ),  cf.  Arnold  [1964]. 

3.  If  g  =  g(x)  is  not  explicitly  time  dependent  then  we  have, 
using  Green's  theorem 


fCq°(t  -  tQ)) 


g(q°(t  -  tQ))dt 


( f  1 92  "  f 29-i  )  dt 


(ggtu^v0)^  -  g1(u°,v°)vu)  dt 


j  g2(u,v)  du  -  g^u.vjdv 
trace  Dg(x)  dx. 

Mnt  T 


(3.14) 


Thus  the  formula  obtained  in  Andronov  et  al  [1971]  is  a  special  (planar) 
case  of  the  more  general  Melnikov  function  which  describes  the  'splitting' 
of  the  perturbed  saddle  seperatrices. 

We  now  turn  to  the  case  in  which  the  pert jrbation  g  ®  g(x,c;u) 
depends  upon  parameters  y  €  IR^ .  For  simplicity  we  take  p  °  1. 


Corollary  3.4.  Consider  the  parametrized  family  x  =  f(x)  ♦  cg{x,t;u), 

u  t  IR  and  let  hypotheses  Al - 3  hole.  Suppose  that  the  Melnikov  function 

2 

H(tA,yl  has  a  quadratic  zero  H(?,^)  =  (r,ub)  B  0  but  ^  (i.y  )  f  0  and 

-  O  it0 

-U 

~  (t,yt)  f  0.  Then  *  0(c)  is  a  bifurcation  value  for  which 

w  — — —  &  Q  '  "  ***  **  ■■ 


quadratic  homed inic  tangencies  occur  in  the  family  og  systems. 
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Proof.  By  hypothesis,  using  (3.5),  we  have 

d(tQ,u)  3  e(a(p  -  ub>  +  S(tQ  -  t)2}  +  0(e[y  -  yb|2)  +  0(e2)  (3.15) 

where  we  have  expanded  in  a  Taylor  series  about  (tQ,y)  *  ( x » )  and  a, 

S  are  finite  constants.  Taking  e  sufficiently  small  we  find  that 
d(t  ,y)  has  a  quadratic  zero  with  respect  to  tQ  for  some  yg  near  ub. 
and  hence  that  Wu(pT),  Ws(pT)  have  a  quadratic  tangency  near  q°(0)  on 

C  £ 

z\  a 

This  result  is  important,  since  it  permits  us  to  verify  in  specific 
examples  one  of  the  hypotheses  of  Newhouse's  [1979]  theorem  on  wild  hyper¬ 
bolic  sets  (see  Section  5). 

We  next  turn  to  the  perturbations  of  the  per'  -‘i c  motions  qa  inside  r. 


3.2  Bifurcations  of  Subharmonics 

Once  more  we  start  ~ith  a  perturbation  lenroa: 

Lerrma  3.5.  Let  qa(t  -  t  )  be  a  periodic  orbit  of  the  unperturbed  system 
t  0 

based  on  z  0  wi th  period  T  .  Then  there  exists  a  perturbed  orbit 
q^(t,tQ),  not  necessarily  periodic,  which  can  be  expressed  as 

•  qa(t  -  tQ)  ♦  tq^(t.t0)  ♦  D(c2)  (3.16) 

uniformly  in  t  €  [tQ,t0  *  T  ],  for  c  sufficiently  small  and  all 

o  €  (-1,0). 

Proof.  The  proof  relies  heavily  upon  the  geometrical  structure  of  the  per¬ 
turbed  stable  and  unstable  manifolds  established  in  lema  3.2.  We  again 
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fix  a  neighborhood  and  take  a  curve  of  initial  conditions  qa(Q)  G  E  0 

not  Ivino  in  LI  and  with  Tim  qa(0)  =  q°(0).  Any  orbit  qa(t  -  tj 
v  a-*  0  0 

starting  on  such  a  curve  takes  a  finite  time  to  reach  the  boundary  of 

as  t  increases  or  decreases  and  hence  we  have  !qa(t,t  )  -  qa(t  -  t  )|  = 

0(c)  for  t  €  (t  -  tp  tQ  +  t2),  say.  Once  in  Uy  the  perturbed  and 

unperturbed  orbits  may  take  arbitrarily  long  to  pass  through  and  exit, 

since,  as  a  -*■  Q,  they  pass  arbitrarily  close  to  the  saddle  point  pQ  (or 

to  y£) •  However,  for  fixed  e  =  eQ  we  can  find  a  set  of  orbits  lying 

sufficiently  close  to  the  stable  and  unstable  manifolds  which  remain  within 

0(e)  of  those  manifolds  until  they  enter  ace  neighborhood  U  of  p, 

G£ 

which  is  chosen  to  contain  y  ,  Figure  3.2.  This  implies  that  we  can 

extend  our  estimate  uniformly  to  t£  (tQ  -  t^,  tQ  +  t^),  where  tj  +  t^ 

is  the  time  required  for  the  unperturbed  orbit  qa  to  pass  from  the 
boundary  of  and  return  to  it.  It  remains  to  check  that  qa(t,tQ) 

and  qa(t  -  t  )  remain  within  0(e)  during  the  arbitrarily  long  passage 


Figure  3.2.  Orbits  i;.  and  UC£ 
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through  U  .  This  follows  from  the  ‘shear*  on  the  flow  near  p;  one  has 

to  check  that  qa  exits  from  U  'at  the  correct  time1  and  since  orbits 
e  ce 

passing  near  y£  take  arbitrarily  long  to  exit,  and  those  near  the  boundary 
of  U  arbitrarily  short  times,  at  least  one  orbit  can  be  found  for  any 

CE 

given  (unperturbed)  time  of  passage.  It  follows  that  an  initial  condition 
t 

q“(t0)€E  0  can  be  picked  with  0(e)  of  qa(G)  and  q°(0)  such  that  the 
orbit  q^(tQ)  will  remain  within  0(e)  of  qa(t  -  tQ)  (and  q*(t  -  tQ) ) 
until  it  reaches  U.  .  It  will  then  'transfer  its  allegiance'  to  qu(t  -  t  ) 
until  it  once  more  reaches  an  e-neighborhood  of  qa(0).  Throughout  it  re¬ 
mains  within  0(e)  of  qa(t  -  t  ).  This  takes  care  of  orbits  with  un¬ 
perturbed  periods  T^  larger  than  some  T  =  T^(eQ)  depending  on  eQ. 

For  orbits  with  periods  shorter  than  V(e0)  a  standard  Gronwall  estimate 

ensures  e-closeness.  Then  we  have  our  result  for  all  T  and  e  *  e„. 

a  o 

But  since  f  and  g  are  Cr  the  solution  will  vary  smoothly  in  e  and 
thus  the  result  holds  for  all  0  <  e  £  eQ.  For  details  see  Greenspan  [1981 ].■ 

ct 

We  next  define  the  subharmonic  Melnikov  function..  Letting  q  (t  -  tQ) 

mT 

be  a  periodic  orbit  of  period  — ,  with  m  and  n  relatively  prime,  we 
set 

mT 

Km/n(t0)  =  j  f(qa(t  -  tQ))  ~  g(qa(t  -t0),t)dt.  (3.17) 

Theorem  3.6.  If  (t  )  has  simple  zeros  and  maxima  and  minima  of  0(1), 

and  dTft/dha  0,  then  for  0  <  e  <_  e(n)  (3.1)  has  a  subharmonic  orbit  of 
period  mT.  Jjf  n  =  1  then  the  result  is  uniformly  valid  in  0  <  e  £  e  . 


Proof.  A  calculation  similar  to  that  of  Theorem  3.3  shows  that 
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,mT 


f(qa(0) )  -  (q®(tn  +  mT,tn)  -  q“(t0,tQ))  *  e  _  f(qa(t 


V> 


*  g{qa(t  -  tQ) ,t)dt/ 1 f (qa(0) ) I  +  0(e2).  (3.18) 


Thus  if  Mm//n(t  )  has  a  simple  zero  then  there  is  a  perturbed  orbit 
0  t 

q“(t,t  )  which  leaves  q^(tQ)  and  returns  to  £  0  at  q“(tQ  +  mT)  with 

the  vector  q“(tQ  +  mT)  -  q“(tQ)  C  Z  0  parallel  to  f(qa(0));  i.e.  the 

projection  onto  fL(qa(0))  vanishes.  Letting  M^(tQ)  =  j  f(q^)  -  g(q^,t)dt 

for  8  near  a,  it  is  clear  that  depends  smoothly  on  8.  For 

8  <a ,Te  <  T  and  for  8  >  a.T.  >  T  ;  it  therefore  follows  that  we  can  find 
'  p  w  0  8  p  a 

perturbed  orbits  q^ ,  q£2,  8^  <  a<  82  such  that  the  vectors 

8.  8.  t  8. 

q  (t  +  mT)  -  q  ](t  )  C  z  0  are  parallel  to  f(q  ’(0)),  but  that  they 

£  v  GO 

have  opposite  orientations.  Thus  the  curve  of  initial  conditions  connecting 

8,  8?  t 

q  '(t0)  to  qe  (tQ)  is  mapped  back  to  the  section  z  0  under  m  iterates 

t  t 

of  P  0  as  indicated  in  Figure  3.3.  It  follows  that  (P,°)m  has  a  fixed 

G  fc» 

point  near  qa(0),  and  hence  that  there  is  a  subharmonic  of  order  m/n. 

The  nonuniformity  for  n  >  1  arises  because  Lemma  3.5  applies  only  to 
orbits  with  times  of  length  T^  =  mT  making  one  pass  through  U^(p)  and 
'ultrasubharmonics '  of  period  mT/n  make  n  passes  through  Uy(p).  ■ 


•fixed  point 


Figure  3.3.  The  proof  of  Theorem  3.6. 
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We  also  have  a  bifurcation  result  analogous  to  Corollary  3.4: 


Corollary  3.7.  Consider  the  parametrized  family  x  =  f(x)  +  eg(x,t;u). 


u  €  IR,  and  let  hypotheses  Al-3  hold.  Suppose  that  Mm/,n(t0,y)  has  a 


quadratic  zero  M 


,m/n  _  3M 


im/n 


3tn  =  °» 


3^  n 

"Vti: 


3M 


m/n 


t  0  at 


V 


Then 


um/n  =  +  a  bifurcation  value  at  which  saddle-nodes  occur. 


Proof.  The  proof  is  identical  to  that  of  Corollary  3.4.  We  will  consider 
the  stability  types  of  the  orbits  created  in  this  bifurcation  in  Section  3.3.H 

The  next  result  is  a  generalization  of  oneobtained  by  Chow,  Hale  and 
Mallet  Paref.  [1980].  It  implies  that  the  homoclinic  bifurcation  is  the 
limit  of  a  countable  sequence  of  subharmonic  saddle- node  bifurcations. 

Theorem  3.8.  Let  Mm//1(tQ)  =Mm(to),  then 

lira  Mm(t  )  =  M(t  }.  (3.19) 

m-*<»  0 

Proof.  We  must  show  that  the  integral 

ml 

m  r? 

M  (tQ)  *  j  f(qa(t  -  t0)5  -  g(ga(t  -  tQ),t)  dt  (3.20) 

_  mT 

T 

converges  to 

<x> 

M(t0)  =  f  f(q°(t  -  t  ))  -  g(q°(t  -  t  ),t)  dt  ,  (3.21) 

*  -CO  ^ 

as  a  -*■  0  and  m  ».  (Note  that  the  periodicity  of  Mm(tQ)  implies  that 
we  can  change  the  limits  from  0  +  mT  to  -  —--.Jy.)  Letting 


ra  =  {qa(t)|t  6  [0,Ta)}  and  r°  =  {q°(t} |t  €  IR}  u  {pQ}  we  select  a 
neighborhood  Uy(p)  such  that  the  arc  length  of  r°,  ra  n  Uy(p)  is  less 
than  v.  Choose  t  such  that  q°(-t)  and  q°(x)  both,  lie  within  U  . 
Then  for  a  close  enough  to  zero,  qa(±x)  also  lie  in  U  .  We  have 


M(t0}  -  rf"(t0)  * 


f(q°)  A  g(q°,t)dt 


f(qa)  a  g(qa,t)dt 

'  -X 


+ 


f  f(q0)  A  gCq°»t)dt  + 


f(q°)  *  g(q°,t)  dt 

'T 


(3.22) 


-  f  f(qa)  ^  g(qa»t)dt  - 
>  mT 

“T 


mT 


rT 

f(qa) 


A  g(qa,t)dt 


The  smoothness  of  f  and  g  and  continuity  of  solutions  with  respect  to 
initial  conditions  implies  that,  given  v  >  0,  there  is  an  a1  <  0  such 
that,  for  a  e  (a',0)  the  first  [bracketed]  term  of  (3.22)  is  less  than 
v.  Clearly,  as  a'->  0”,  m  00 .  The  second  term  may  be  expressed  as  the 
difference  between  two  integrals  over  arcs  of  r°  and  Ta.  Using  the  arc- 

length  increment  ds  =  +  v0^  dt  =  |f(q°)jdt  on  F°  and  ds  =  |f(qa)|dt 

on  ra,  the  second  term  becomes 


f(q°)  -  g ( q° , t ) 


ds 

I  f  (q°)  I 


fqa(x) 

c  f(qa)  -  g(qa,t) 
qa(-x) 


ds 

I  f  (qa)  I 


(3.23) 


Our  assumptions  on  g  imply  that  sup  |g(qa,t)|  *  K  <  <»,  and  thus 

qe  Bv(p) 

t  €IR 

that  the  second  term  is  bounded  by  2Kv .  Hence,  for  a  G  (a‘,0) 


3.14 


|M(t0)  -  ^itQ)\  <  (2K+  1)v 
and  |M(tQ)  -  Mm{t0)|  -*-0  as  ot  -*■  0.  ■ 


3.3  Higher  Order  Terms  and  Stability 

We  next  develop  a  perturbation  method  which  enables  us  to  compute  the 

t 

global  structure  of  the  perturbed  Poincar6  map  P£0,  and  to  determine  how 
the  sets  of  subharmonics  and  homoclinic  orbits  are  related.  Our  starting 
point  is  Melnikov  [1963,  §7],  although  we  have  somewhat  modified  his  trans¬ 
formations. 

Since  the  unperturbed  system  is  Hamiltonian,  a  symplectic  change  of 
coordinates  to  action  angle  variables  can  be  found  in  the  interior  of  r°: 


I  =  I(U,v)  0  *  0(U,v) 

Under  this  change  of  coordinates  (3.1)  becomes 


(3.24) 


I 


fai  „  .31 

3u  91  3V  g2 


def 


£F(I,0,t) 


(3.25) 


0  =  Q(I)  +  e 


\ 30  _  ,  30  . 

3u  91  3v  92 


^  G(I,0,t) 


where  H(Ia)  =  tj-  (Ia)  =  v-  is  the  angular  frequency  of  the  unperturbed 

3  a 

orbit  qa(t  -  t  )  with  action  Ia  =  I(qa).  We  now  consider  small  perturba- 

rr\T 

tions  of  a  resonant  orbit  T  =  —  .  Letting 

an  3 


I  =  Ia  +  v£h 

0  =  ft(Ia)t  +  <p  =  t  +  $  ^  ftat  +  <(>,  (3.26) 

we  obtain 


h  s  v£F(Ia,flat  +  <j>,t)  +  eF'(I  »fl  t  +  <j>,t)h  3/2 

+  0[e'^)  , 

i  =  ^ai'd^h  +  e(GCla,na  t  +  *,t)  +  n"Cla)h2) 


3.15 


(3.27) 


where  '  denotes  gj.  Here  we  have  expanded  in  Taylor  series  and  used  the 
fact  that  fl*  7s  0,  since  dT^/dh^  f  0.  Noting  that 


3I.3TaH  =  .  1  f  81  -  1  f 

3u  3H  3u  0(17  f2  and  3v  o[TT  fl 


the  leading  term  of  (3.27)  can  be  rewritten 


h  =  /e  -j-  f(qa(t  -  <f>/Qa) )  -  g(qa(t  -  4>/fla),t) 
fl 

*  =  /efl'(Ia)h. 


(3.28) 


Provided  that  fl'(Ia)  is  bounded,  for  v'e’  sufficiently  small  the  averaging 
tneorem  (cf.  Hale  [1963])  can  be  applied  to  (3.28)  to  yield 


h 


„a 

fl 


1 

mT 


wnT 

•0 


f(qa(H>/fla)) 


a  g(qa(t  -  4>/fla)  ,t)dt 


or 


h.  ^N”/n 


_1 

rS^ 


=  /£  fl'  (Ia)F  . 


(3.29) 


Under  the  averaging  theorem,  the  hyperbolic  or  elliptic  fixed  points  of 
(3.29)  correspond  to  small  periodic  motions  of  (3.27)  and  hence  to  sub¬ 
harmonics  of  order  m/n  of  (3.1).  It  is,  of  course,  no  coincidence  that 
necessary  and  sufficient  condition  for  the  existence  of  such  fixed  points 
is  ’’hat  the  Melnikov  function  Mm/,n  have  simple  zeros  and  that  fl'(Ia), 

fl(Ia)  f  0  (dT  /dh  t  0). 
a  a 

We  note  that  (3.29)  is  a  structural ly  unstable  Hamiltonian  system  with 


3.16 


Hamiltonian 


H 


(3.30) 


where  V(<j>)  3  j Mmy/n(^7fla)d(j),  and  thus  to  determine  the  stability  and 
the  global  behavior  of  orbits  of  the  unperturbed  system  near  the  resonant 
orbit  qa,  we  must  investigate  the  terms  of  0(e).  Letting  f  -  g  «  f*m^n 
+  F(.<j»,t)  where  F  has  period  T  and  zero  mean,  the  averaging  transformation 
is 


h  h  +  *€ 


F(4>,t)dt;  $  $  , 


(3.31) 


where  the  antiderivative  is  defined  up  to  a  t-i ndependent  term,  generally 
taken  to  be  zero.  Using  (3.31),  (3.27)  becomes 


F  3  Ji  ^fMm/n(4./fia)  +  e((F,(Ia,Oat  +  <ji.t)h  -  ^  j 


F  dta'h), 


4)  =  Viii'h  +  e(0"hc  +  G(Ia,aat  +  4>,t)  +  ft*  F  dt)  . 


+  0(e3/2) 


(3.32) 


Since  F  has  zero  mean  (it  is  simply  a  sum  of  fburier  components),  Jf  and 
■—  J  F  also  have  zero  mean  and  on  a  second  application  of  averaging  to  the 
0(e  )  terms  of  (3.32)  we  obtain  (dropping  the  bars) 


h  =  /£  ■^Mm'n(4>/«a)  +  eF'U)h, 
<p  =  4-fl'h  +  e(a“’h2  +  G(4>) )  , 


+  0(e3/2)  , 


(3.33) 


*  where  P",  G  are  the  averages  of  F‘  and  G.  As  Horosov  [1973]  notes, 
this  second  order  averaging  generally  suffices  to  determine  the  stability 
of  the  fixed  points  and  hence  of  the  bifurcating  subharmonics,  at  least  for 
a*  <  ®  and  c  sufficiently  small.  However,  as  we  shall  see  in  our  appli¬ 
cation  in  sections  4  and  6,  one  can  also  obtain  global  information  on  the 


Poincard  map  by  considering  the  time  T  flow  maps  of  the  averaged  systems 
(3.33)  in  the  neighborhood  of  each  resonant  and  nonresonant  periodic  orbit. 
These  results  on  the  full  PoincarS  map  P  follow  from  application  of  the 
averaging  theorem  (Hale  [1969]). 


4.1 


4.  GLOBAL  BIFURCATIONS  OF  PUFFING'S  EQUATION 

In  this  section  we  apply  the  results  of  Section  3  to  the  Duffing  equa¬ 
tion  with  weak  sinusoidal  forcing  and  damping.  Written  as  a  system,  we 
have 

°sv  (4.1) 

v  =  u  -  u  +  e(y  cos  wt  -  6v)  , 

where  the  force  amplitude,  y,  frequency  u>,  and  the  damping  <$  are 
variable  parameters  and  e  is  a  small  scaling  parameter.  For  e  =  0  the 
system  has  centers  at  (u,v)  =  (±1,0)  and  a  hyperbolic  saddle  at  (0,0). 

The  level  set 

2  2  4 

(u,v)  =  -2--^-  +  -^rs0  (4.2) 

is  composed  of  two  homoclinic  orbits, r°,  r°  and  the  point  pQ  =  (0,0).  The 
unperturbed  homoclinic  orbits  are  given  by 

qj(t  -  tQ)  3  (/T  sech(t  -  tQ),-  f£  sech(t  -  tQ)  tanh(t  -  tQ)), 

3  -qj  .  (4.3) 

Within  each  of  the  loops  T°  there  is  a  one  parameter  family  of  periodic 
orbits  which  may  be  written 


where  sn,  cn  and  dn  are  the  Jacobi  elliptic  functions  and  k  is  the 

elliptic  modulus.  As  k  -  1  q^  -  q°  u  (0,0)  and  as  k  *  0  q*  ~  (±1,0).  We 

nave  selected  intial  conditions  at  t  =  t„: 

o 


4.2 


4.3 


M(  tQ  ,y  »(S  ,to) 


00 

/ 

v°(t  -  t  )[y  cos  wt  -  6v°(t  -  t  )]  dt 

-00  0  0 


<x> 

t 


sech(t  -  tQ)  tanh(t  -  t  )  cos  ut  dt 


sech2(t  -  t  )  tanh2(t  -  t.)  dt  . 


(4.10) 


The  integrals  are  easily  evaluated  [the  first  by  the  method  of  residues)  to 
yield 


M(tQ,Y»6.»)  3 


46 

3 


sin  ut  . 


(4.11) 


If  we  define 


R°(w) 


4  cosh 
3i/2t?u 


itu 

T 


(4.12) 


then  it  follows  from  Theorem  3.3  that  if  y/6  >  R°(u),Ws(Pe)  intersects 
WU(P  )  for  e  sufficiently  small  and  if  y/6  <  R°(u),Ws(P  )  nwu(P  )  =  3. 
Moreover,  since  M(to,y,6,u)  has  quadratic  zeros  when  y /6  «*  R°(u),  it 
follows  from  Corollary  3.4  that  there  is  a  bifurcation  curve  in  the  y,  6 
plane  foreach  fixed  u  tangent  to  y  s  R°(u)6  at  y  »  6  ®  0  on  which  quad¬ 
ratic  homoclinic  tangencies  occur.  We  display  some  calculated  bifurcation 
curves  in  figure  4.1,  below. 

We  next  compute  the  Melnikov  function  for  the  resonant  periodic  orbits. 
We  will  only  consider  those  within  T°,  q^{t  -  tQ).  The  resonance  condition 
is,  from  (4.3) 

2K(k) -  k2  »  ~  (4.13) 

and  for  each  choice  of  m,  n  with  (4.13)  can  be  solved  to  give 

a  unique  resonant  orbit  Computing 


4.4 


tf"/n(t „,y,M)  3  I  Vk(n,’n)(t  -  t  )Cy  cos  ut  -  6vk(m‘n)(t  -  tn)]  dt.  (4.14) 


Using  the  Fourier  series  for  (4.4)  and  Remark  3  following  Theorem  3.3,  we 
obtain 

jin/  n  i 


M" 


(tQ ,y,6,uj)  »  -<5Jj(m,n)  -  Y^ro.n.w)  sinu)t0>  (4.15) 


where  U^(m,n)  =  |[(2  -  k2(tn,n))2E(k(m,n))  -  4k,2(m,n)K(k(ra,n))]/(2  -k2(ra,n)) 
and 


J2(m,n,oj) 


j  0  ;  n  t  1 


Here  E(k)  is  the  complete  elliptic  integral  of  the  second  kind  and  k' 

2  2 

is  the  complementary  elliptic  modulus  k‘  =  1  -  k  .  Defining 


J1  (m,l ) 

3  ^TiT^T  ) 


(4.16) 


we  conclude  from  Theorem  3.6  and  Corollary  3.7  that  if  y/6  >  Rm(«)  then 
there  is  a  pair  of  subharmonics  of  order  m  (period  2rrm/ui)  which  appear 
on  a  bifurcation  curve  tangent  to  y  °  Rm(w)6  at  y  M  *  0, 

Routine  computations  verify  that 


lim  ,y,4,w)  =  H(t  ,y,6,w)  , 


(4.17) 


that  the  limit  is  approached  from  below  and  that  the  rate  of  convergence  is 
extiemely  rapid  (cf.  Figure  6.B). 

Similar  computations  for  the  orbits  lying  outside  T°  U  ((0.0)}  U  '° 


give  the  Melnikov  function 


4.5 


^n(VY,6,*)  3  -6J-J  Ctn.n)  -  YJ2(m,n,u>)  sinujt0  (4.18) 

with 

J-j (m,n)  =  |  [(2k2(m,n)  -  1  )4E(k(m,n))  +  4k1 (ro,n)2K(k(m,n)]/(2k2(m,n)  -  1)3/2 
and 


J^ra.n.u) 


0  ;  n  f  1 ,  m  even, 

^  ifh u>  sech  [ ( ra ^ T) ^ ' )  *  n  3  1  >  m  odd. 


In  this  case  we  obtain  a  sequence  of  bifurcation  curves 

a  a  def  /srn  o 

y  »  vh  (mj  ,u>)  R  (o>)6  accumulating  on  y  =  R  (u))<$  from  above. 

Here  k(m,n)  is  the  unique  solution  of  the  resonance  relation 

Tk  -  4K(k)/2k2^  1  =  ^  .  (4-19) 


Formore  details  on  the  above  computations  see  Greenspan  [1981]. 

In  Figure  4.1  we  show  some  of  the  approximate  bifurcation  curves 
y/<5  .  R° (uj ) ,  Rm{w)  for  orbits  both  inside  and  outside  the  separatrix.  Below, 
in  Figure  6.8  we  show  how  the  subharmonic  bifurcation  curves  accumulate  on 
‘.he  homed, nic  uifurcation  curve  as  m  -  «  for  fixed  u>.  We  leave  the 
discussion  of  the  detailed  structure  of  the  Poincar4  maps  for  various  param¬ 
eter  values  (y.o.w)  to  Section  6. 


V 

z> 


S 


Figure  4*1  (a)  Bifurcation  curves  Rm(o))  for  subharmonics  and 
homoclinic  orbits  inside  rQ  +  U  {0,0}  U  VQ  . 
R1  ■  subharmonic  of  order  i,  R^  *  horaoclinic 


[•  J  m  •] 


t. 


7J 


Figure  4«1  (b)  Bifurcation  curves  Rm(w)  for  subharmonics  and 
homoclinic  orbits  outside  Tq  +  U  {0,0}  U  Fq~. 


5.1 


5.  SMALE  HORSESHOES,  NEWHOUSE  SINKS  AND  CHAOTIC  MOTIONS 

In  this  section  we  provide  a  brief  review  of  some  useful  results  in  the 
abstract  theory  of  dynamical  systems,  and  in  particular  of  two  dimensional 
diffeomorphisms.  For  background  material  see  Smale  [1963,1967],  Nitecki 
[1971],  Moser  [1973]  or  Newhouse  [1980]. 


5.1  The  Smale  Horseshoe 


Here  we  outline  an  important  example  of  a  planar  diffeomorphism  due 

to  Smale  [1963].  We  shall  adapt  the  map  to  our  specific  application. 

2 

Consider  a  Hk.p  F  :  Q  IR  defined  on  the  square  Q  »  [0,1]  x  [0,1] 
such  that  F  is  linear  on  the  two  horizontal  strips  =  [0,1]  x  [0,a], 
Ho  =  [0,1]  x  [1  -a, a] 


DF(p)  = 


with 

0  ' 

-6  0  • 

;  p  €  H.  and  DF(p)  = 

s 

1/a. 

l 

Lo  -1/aJ 

;  P€  Ho,  (5.1) 


where  0  <  6  <  a  <  1/2.  F  is  chosen  such  that  the  image  of  the  strip 
[0,1]  x  (a,  1  -  a)  falls  outside  Q  but  F(H^)  —  and  F(H2)  V2 
are  two  disjoint  vertical  strips  lying  in  Q  as  shown  in  Figure  5.1.  By 
our  construction,  F  has  a  hyperbolic  saddle  point  (0,0)  and  F/Q  has 
constant  Jacobian  g/a  <  1  and  in  thus  a  dissipative  map  (Newhouse  [1979]). 
The  nonwanderinq  set  A  c  Q  of  F  is  obtained  by  intersecting  all  forward 
and  backward  iterates: 

A  =  Av  n  V  Av  *  °  \  =  n  F"n(Q)  •  (5-2) 

n>0  nH) 

Thus  =  Q  n  F(Q)  =  V|  u  is  a  pair  of  vertical  strips  each  of  width 
6.  It  is  not  hard  to  see  that  =  h  Fn(Q)  is  2^  disjoint  vertical 

v  n^o 


5.2 


F([0,1]  x  (a,l-o)) 
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/ 

V1 

r 
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Figure  5.1.  The  Horseshoe 


strips  each  of  width  0  and  thus  that  A  =  C  x  [0,1]  is  a  Cantor  set  of 


vertical  intervals.  Similarly  F"n(Q)  is  2k  disjoint  horizontal 

strips  each  of  width  ak  and  A.  =  [0,1]  »C  is  a  Cantor  set  of  horizon- 

K  Ct 

tal  intervals.  Thus  A  =  C0  x  C  is  also  a  Cantor  set. 

3  a 

To  describe  the  orbits  of  F|^  we  assign  to  each  point  x  €  a  a 
bi-infinite  symbol  sequence  {a ^ (x) }  =  {. . .a^a^ *a0a]a2  ••••*'  chosen  such 
that  a.(x)  =  1  (resp.  2)  if  FJ(x)  e  H.  (resp.  ti9)  and  a  .(x)  *  1  (resp. 
2)  if  F'^"^(x)  €  Vi  (resp.  V^).  The  action  of  F  on  A  then  corresponds 
to  the  shift  a  on  the  space  E  of  all  such  symbol  sequences,  and  in  fact 
there  is  a  homeomorphism  $  such  that  the  following  diagram  commutes: 


A 

w 


E  • 


5.3 


...  i  '•■.v;,  ,-s  >■  -ix  „  gWfr*  j  _ _ _ _ _ _ 


To  every  periodic  sequence  there  corresponds  a  periodic  orbit  of  F|A;  for 
example  the  fixed  point  (0,0)  has  a  symbol  sequence  {..  111*111...} 
and  there  is  a  second  fixed  point  with  sequence  {...  222*222  ...}.  For 
details  see  Moser  [1973].  This  method  of  symbolic  dynamics  enables  us 
to  prove  the  following. 

Proposition  5.1.  The  map  F  has  a  countable  set  of  periodic  points,  with 
points  of  arbitrarily  high  period  and  an  uncountable  set  of. non  periodic 
recurrent  motions.  The  periodic  points  are  dense  in  A  and  there  is  an 
orbit  dense  in  A.  All  the  periodic  points  are  of  saddle  type. 

The  last  assertion  follows  directly  from  the  form  of  the  linearized 
map  OF.  it  is  more  difficult  to  prove  that  F j ^  is  structurally  stable, 
but  this  latter  important  fact  implies  that  any  sufficiently  C^-close  map 
F  also  possesses  an  invariant  cantor  set  A  homeomorphic  to  A.  Thus 
the  complicated  dynamics  of  F  cannot  be  removed  by  small  perturbations. 

The  'piecewise  linear'  map  F  represents  an  idealization  of  a  map 
which  is  embedded  in  the  dynamics  of  any  system  with  transverse  homoclinic 
orbits.  This  result,  the  Smale-Birkhoff  homoclinic  theorem  (cf.  Birkhoff 
L 1 9273  Smale  [1963],  Moser  [1973])  follows  from  the  lambda  lemma  and  the 
homomorphic  structure  of  Figure  2.2,  although  it  can  be  proved  independently. 
We  select  a  small  rectangle  S  bounded  by  pieces  of  the  stable  and  unstable 
manifolds  W  (p),  W  (p),  with  a  transverse  homoclinic  point  q  at  one 

"*i 

corner.  There  are  then  integers  £, ,  £?,  n.  such  that  A  =  P  (S), 
h  i  a  i 

8  =  P  t$)  are  disjoint  and  lie  in  any  neighborhood  Uv(p),  and  the  map 

ni  ni 

P  :  B  +  A  is  well  approximated  by  its  linearization  DP  (p).  This  en¬ 
ables  one  to  find  horizontal  strips  c  B  whose  images  under 


JUH+n.  (JL+JL+n.)  •(i1+S,?+n?) 

P  L  1  are  vertical  strips  V-jjVgC  ^he  maP  ^ 

then  has  a  horseshoe.  A  careful  study  of  the  linearized  map  enables  one  to 

obtain  estimates  of  hyperbolicity  and  show  that  the  behavior  of 
( £«  +  2,+n. )  ( £.. +2.p+n . ) 

p  4  1  is  homeomorphic  to  that- of  the  'idealized'  map  F. 

Specifically,  for  our  problem  we  have 


t  t  t 

Theorem  5.2.  If  the  Poincare  map  P£°  :  E  0  -*•  E  0  possesses  a  transverse 


homoclinic  point  q  to  a  hvDerbolic  saddle  point 


then,  in  a  neiah- 


borhood  of  q  ,  some  power  (P  °)M  possess  anjinvariant_zero-dimensional 

£  +  t 

cantor  set  A  on  which  (P  °)  is  conjugate  to  a  shift  on  two  symbols. 


Corollary  5.3.  (P  0 )  /A  possess  a  dense  set  of  periodic  points,  there  are 
pints  of  arbitrarily  high  period  and  there  is  a  non-periodic  orbit  dense 

A 

in  A.  Moreover  A  is  structurally  stable  in  the  sense  that  if  P  is 

t 

0  A  M 

sufficiently  C'  close  to  P_ ,  then  (P)  has  a  non-wanderinq  set  A 

e  * 

^  n  -N 

and  there  is  a  homeomorphism  h  such  that  (P  )  °h  *  L°P  . 


The  horeseshoe's  non-wandering  set  A  is  extremely  complicated.  It 
is  of  saddle  type  and  has  one  dimensional  stable  and  unstable  bundles 
WS(A),  • WU(A)  which  are  locally  the  products  of  Cantor  sets  and  curves. 
Letting  Ws(Q,0),  Wu(0,0)  be  the  stable  and  unstable  manifolds  of  the 
fixed  point  {...111*111...}  of  the  idealized  horseshoe,  F,  we  have 
WS(A)  3  C2,(WS(0,0))  and  WU(A)  =  Ci(Wu(0,0) ) ,  and  since  such  a  A  can 
be  chosen  arbitrarily  close  to  any  transverse  homoclinic  point  q£  and 
hen'-e  to  the  saddle  point  p£  in  our  application,  it  follows  that 

WS(A)  =  Ct(Ws(p£)) 

WU(A)  =  Cz(Wu(p£)). 


«La> .c.-  .  MttHt—a*.  —  nm-AWfc. 


Thus  the  structure  of  the  stable  and  unstabl e  manifolds  Ws(p  ),  Wu(p  )  Is 

C  £ 

doubly  Important  not  just  for  establishing  the  existence  of  horseshoes,  but 
also  in  their  structure.  We  shall  return  to  this  in  Section  6. 

5.2  Homoclinic  Tanqencies  and  Newhouse  Sinks 

We  now  turn  to  the  bifurcations  in  which  horseshoes  are  created.  Newhouse 
[1974,  1979,  1980]  has  studied  the  situation  and  has  shown  that  the  lack  of 
hyperbolicity  which  occurs  when  stable  and  unstable  manifolds  havea  quadratic 
tangency  persists  under  small  perturbations.  We  first  require  a  definition. 

A  hyperbolic  (basic)  set  A(F)  for  a  two  dimensional  diffeomorphism  F  is 
called  wild  if  there  is  a  C  -neighborhood  ll  of  f  such  that  for  any  6  6N 
there  are  points  x,y  €  A( G)  such  that  Wu(x)  and  Ws(y)  are  tangent  some¬ 
where. 


Theorem  5.4.  (Newhouse  [1979]).  Let  p  be  a  dissipative  saddle  point  of 

a  two  dimensional  Cr  diffeomorphisr,.  F  (|det(DF(p))|  <  1).  Suppose 

us  r 

W  (p),  W  (p)  are  tangent  at  some  point  q.  Then  arbitrarily  C  near  F 

there  is  a  diffeomorphism  G  naving  a  wild  hyperbolic  set  near  the  orbit  of 

q. 

Wild  hyperbolic  sets  are  important  because  they  imply  the  existence  of 
countably  many  stable  periodic  orbits. 

Tneorem  5.5.  Suppose  F  is  a  two  dimensional  Cr  diffeomorphism  with  a 
wild  hyperbolic  set  a  containing  a  dissipative  periodic  point.  Then  there 
is  a  neighborhood  w  of  F  and  a  residual  subset  8  C  N  such  that  if 
G  C  8.  then  G  has  countable  many  periodic  sinks. 
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For  a  proof  of  Theorem  5.5  see  Newhouse  [1980].  One  step  in  the  proof 
closely  resembles  a  result  obtained  independently  by  Gavrilov  and  Gilnikov 
[1972,1973],  who  showed  that,  for  certain  two  dimensional  di f feomorphi sms , 
the  homoclinic  tangencies  were  the  limit  of  an  infinite  sequence  of  saddle- 
node  bifurcations  in  which  pairs  of  periodic  sinks  and  saddles  are  created. 

To  illustrate  their  argument,  consider  a  one  parameter  family  of  maps 

such  that  for  y  <  0  Wu(p)  nws(p)  =  0  while  for  y  >  0  we  have  transversal 

intersections  and  at  y  =  0  FQ  has  a  quadratic  tangency  as  in  Theorem  5.4. 

We  pick  a  coordinate  system  (x,y)  such  that,  in  a  neighborhood  Uv(p), 

Ws(p)  and  Wu(p)  are  given  by  y  =  0  and  x  *  0  respectively.  Let 
q^  =  (0,y)  and  =  (7,0)  =  F^q^)  be  (non- transverse) homoclinic  points 
and  Ap  A,,  be  neighborhoods  of  q-| ,  lying  in  as  shown  in  Figure 
5.2.  The  map  Fk  :  A^  -*  is  well  approximated  by 

fXK  01 

DF(p))k  .  k  ,  (5.3) 

LC  yK  j 

where  0  <  \  <  1  <  y  and  \y  <  1  and  k  ■*  m  as  we  pick  domains  in  A^ 
closer  and  closer  to  Ws(p).  The  form  of  the  map  f£  :A2  ■"  Ai  (n  *s  fixed) 
may  also  be  approximated  since,  by  hypothesis,  the  vertical  line  Wu(p)  nA 

2 

is  mapped  to  the  parabola  Wu(p)  n  A-j  and  the  parabola  Ws(p)  n  A2  to  the 
horizontal  line  W  (p)  n  A^ .  When  y  j*  0  the  images  also  undergo  a  verti¬ 
cal  translation.  Using  these  facts,  and  normalizing  the  maps  to  remove  re¬ 
dundant  parameters,  we  arrive  at  an  approximate  quadratic  map 

Fu  k  *  fy,k  :  (x  "*•  y>  *  (*^ky  -  7),  Akx  +a(Yky  -  y)2  -  u)  .  (5.4) 

(We  have  selected  the  orientation  such  that  a  >  0).  It  is  interesting  to 
note  that  this  is  essentially  the  same  quadratic  map  studied  by  Henon  [1976] 
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Figure  5.2.  Quadratic  homoclinic  tangencies 


and  conjectured,  on  the  basis  of  numerical  evidence,  to  possess  a  strange 
attractor. 

We  fix  a  horizontal  strip  H(k)  C  A1  and  an  integer  k  such  that 
F  (H(k))  rt  Ag  is  a  vertical  strip  and  hence  Fn+l4(K(k))  n  A-j  is  non  empty. 
Equation  (5.4)  then  permits  us  to  compute  the  approximate  bifurcation  value 
for  which  a  pair  of  periodic  orbits  of  period  n  ♦  k  appear  in  a  saddle  node 
bifurcation: 

„  ,k-  -k-  1  ,,k  -k.2  c. 

uk  =  Xx-y  y  -  U  ♦  Y  )  •  (5.5) 


It  is  easy  to  check  that  uk  -*0*  as  k  -  **,  showing  that  the  homoclinic 
tangency  at  u  °  0  is  the  limit  of  a  countable  sequence  of  saddle  node 
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bifurcations.  For  sufficiently  large  k,  these  subhartnonics  of  order 
m  n  +  k  are  precisely  those  found  by  Melnikov's  method  in  Sections  3 
and  4,  since  the  perturbation  methods  used  there  yield  all  subharmonics 
which  pass  once  through  U  . 

The  two  fixed  points  of  (5.4)  existing  for  y  >  y^  are  a  saddle 
and  a  sink  and  another  straightforward  analysis  shows  that  the  sink  bifur¬ 
cates  to  a  sink  of  period  two  at 

^  =  Xk7  -  y*ky  +  (Ak  +  y“k)2  ,  (5.6) 

leaving  to  an  orbit  of  period  2m  for  the  map.  The  bifurcation  sequence 
y£  also  accumulates  on  y  =  0  from  below.  Recent  work  of  Lanford  [1981] 
has  shown  that  quadratic  maps  of  the  plane  with  small  Jacobian  determinant 
undergo  countable  sequences  of  period  doubling  bifurcations  precisely  as  do 
one  dimensional  maps  to  which  they  are  sufficiently  close  (Feigenbaum  [1978], 
Collet,  Eckmann  and  Lanford  [1980]).  The  present  map  has  } Dg  .  |  a  (xy) 

V  f  K 

and  therefore  falls  into  this  class  for  k  large.  We  conclude  that  y£ 
is  but  the  first  of  a  countable  set  of  period  doubling  bifurcations  which 
must  take  place  before  y  3  0,  since  for  y  >  0  there  is  a  full  horseshoe 
in  Aj  near  q^,  containing  points  of  all  periods  for  the  map  F  .  In 
fact  countable  sets  of  further  bifurcations  must  also  occur,  in  which  points 
of  period  2J-3(n  +  k),  2^-5(n  ♦  k);  j  *  1,  2,  ...  etc.  are  created.  Since 
such  sequences  occur  for  each  sufficiently  large  k  we  have  a  countable 
sequence  of  countable  sequences  of  bifurcations  all  accumulating  on  y  a  0. 
One  can  thus  find  sinks  of  arbitrarily  high  period  in  any  neighborhood  of  a 
non  transverse  homoclinic  point  q,  and  finite  sets  of  these  sinks  persist 
when  the  homoclinic  tangency  is  broken.  Newhouse's  results  [1979,1980]  give 
more  irforraation  and  in  particular  show  that,  in  the  case  of  tangencies,  the 
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closure  of  the  set  of  homoclini.c  points  of  a  suitably  perturbed  map  is  con¬ 
tained  in  the  closure  of  the  set  of  periodic  sinks.  This  implies  that  any 
homoclinic  orbit  can  be  approximated  arbitrarily  closely  (for  finite,  but 
arbitrarily  long  times)  by  a  stable  periodic  motion  and  thus  the  horseshoe 
might  become  'physically  observable'.  On  the  basis  of  this  work  Newhouse 
[1979]  conjectured  that  some  of  the  numerically  observed  strange  attractors 
(such  as  that  of  Hdnon  [1976])  might  in  fact  be  orbits  lying  in  the  stable 
manifolds  of  long  period  sinks.  The  question  is  still  open,  and  we  shall 
return  to  it  in  the  next  section. 

Before  closing  we  note  that,  even  without  Theorem  5.4,  the  existence 
of  a  single  parameter  value  y  at  which  a  homoclinic  tangency  occurs  im¬ 
mediately  implies  a  countable  set  of  such  values,  since  the  stable  and  un¬ 
stable  manifolds  accumulate  on  one  another,  by  the  lambda  lemma.  Thus  the 
'first'  tangency  is  immediatley  followed  by  infinitely  many  others,  each 
with  their  attendant  saddle-node  and  period  doubling  bifurcations. 
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6.  THE  GLOBAL  STRUCTURE  OF  SOLUTIONS  OF  PUFFING'S  EQUATION 

In  this  final  section  we  make  use  of  the  bifurcation  and  perturbation 
computations  of  Sections  3  and  4  and  the  abstract  theory  of  Section  5 to  obtain 
a  partial  understanding  of  Duffing's  equation  (4.1)  for  £  small.  We  start 
by  considering  the  phase  portraits  of  the  averaged  systems  (3.29)  and  (3.33) 
near  each  resonant  periodic  orbit  inside  the  separatrix  r°.  Similar  conclu¬ 
sions  apply  within  r°  and  outside  the  appendix. 


6.1  Subharmonics  and  Their  Domains  of  Attraction 

The  material  in  this  and  the  following  section  is  deduced  directly  from 
the  theory  of  Section  3,  the  results  of  Section  4,  and  standard  results  in 
dynamical  systems  theory  such  as  that  reviewed  in  Section  5.  In  these  two 
sections  we  collect  various  proven  facts  oonceming  Duffing's  equation,  al¬ 
though  we  do  not  state  precise  theorems.  In  the  concluding  Section  6.3  we 
go  on  to  make  some  conjectures  which  we  feel  have  a  firm  basis  in  these  facts. 
Within  r°  the  action  transformation  gives 

Ik  -  I(k)  »  ~  [2(2  -  k2)E(k)  -  4k'2K(k}]/(2  -  k2)3/2  ,  (6.1) 

where  the  elliptic  modulus  k  =  depends  implicitly  upon  u  and  v  through 

v2  2  4  k2-l 

tne  relationships  H(u,v)  =  -n-  -  ♦  V  s  - ~T~7  -  Using  (4.13)  we  have 

1  2  4  (2  -  k V 

r''(,k>‘3f  (tttt);  k’-SsT(kl/3rI(k> 

‘i-i  k 


.  »Z(2  -  kZ)f(2  - /)E(k)  -  U  >  fc2)f2K(k)J  . 
4k  VW 


(6-2) 


b.Z 


where  k  3  k(m)  is  defined  by  the  relationship  K(k )Jl-  k2  =  um/u  on  a 

2  2 

resonant  orbit  of  period  Imj m  and  k1  3  1  -  k  .  Thus  (3.29)  and  Its 

integral  (3.30)  become 

<J>  3  ^'(mjh  ,  (6.3) 

h  3  Ji  (-<S0-|(m)  -  ^(m)  si n(m4>) )  , 


and  H  3  Jz 


Q!M  h2  +  J_ 

2  n  2tt 


6J.j  (m)6<J)  -  y 


Oo(m,w) 


m 


cos  (mcp ) 


(6.4) 


where  we  have  used  equation  (4.15)  (with  J^(m,l)  —  J^(m))  and  the  fact 
that  Q{la)  =  w/m  at  resonance.  Noting  tnat  6J.(m),  y^111*'13)  are  POs*t-ive 
while  n‘(m)  is  negative  we  easily  sketch  the  phase  portraits  of  (6.3)  for 
the  three  cases  y/6  <,  3  ,>  Rm(w),  as  in  Figure  6.1.  Note  that,  for 
y/6  >  Rm(u>)t  in  each  resonant  band  there  are  m  saddles  and  m  ‘centers' 
corresponding  to  a  pair  of  subharmonics  of  order  m  in  the  original  system. 

To  study  the  actual  stability  type  of  the  fixed  points  which  appear  as 
centers  in  the  O(v^)  approximation,  we  compute  the  averaged  higher  order  terns 
F‘,  6  of  Eq.  (3.33).  After  tedious  but  routine  computations  we  obtain 


ptxJC'tt 

where  and  K,  are^conscants  found  by  Integrating  products  of  elliptic  and 
trignocetric  function?,  much  as  in  the  first  order  terras  computed  in  section 
four.  We  noce  that  we  do  not  need  Xi  explicitly,  and  that  tne  trace  (up  to 
0(e))  is  given  by 


<cf .  eg .  3. 33)) . 


These  fixed  points  are  thus  found  to  be  sinks,  and  the  bifurcations  in  which 
they  are  created  simple  saddle  nodes,  rather  than  the  doubly  degenerate 
Hamiltonian  bifurcations  occuring  in  (6.3)  when  y/6  =  Rm(ui).  We  show  the 
phase  portraits  in  Figure  6.2.  The  stable  manifolds  or  domains  of  attrac¬ 
tion  of  the  sinks  are  shaded. 

The  methods  used  here  are  quite  standard  and  include  a  check  on  the 
non  existence  of  closed  orbits  by  Bendixon's  criterion.  Morosov  [1973], 
[1976]  obtained  similar  results  for  a  related  equation. 
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Figure  6.2.  Phase  portraits  of  Eq.  6.3  with  0(e)  terms  added,  m  =  2. 


It  is  important  to  realize,  however,  that  unlike  the  existence  Theorem  3.6 
these  results  are  not  uniformly  valid  since  the  factor  ft'(m)  in  (6.3)  be¬ 
comes  unbounded  as  m  *  in  fact  using  the  asymptotic  behavior 
K(k)  -  ?.n(4/k‘)  as  k  l(k'  -*■  0),  we  have 


a‘(m) 


2 

_ -IT 

4k,2[jin(4/k,)i3 


fe3  e27rm/oj 


irm 


(6.6) 


as  m  -*■  ®.  Thus  averaging  is  valid  for  successively  smaller  regions 
0  <  e  £  e0(mj  as  m  increases,  and,  while  the  invariant  manifolds  and 
periodic  points  of  the  Poincard  map  P^.  are  qualitatively  similar  to  those 
of  Figure  6.2  in  any  given  resonance  band  ,  we  cannot  carry  our  stability 
results  uniformly  to  the  limit  of  homoclinic  bifurcations  as  in  Theorem  3.8. 
Our  results  are,  however,  good  for  subharmonics  of  all  finite  orders  m  <  <*> 
and  hold  within  each  resonance  band  near  the  unperturbed  resonant  orbit  of 
period  2Ttm/u). 

Noting  that  J^m.n)  >  0  while  =0  for  n  f  1  (Eq.  4.15) 

we  conclude  that,  betvieen  the  resonance  bands  all  orbits  decay  ’inward' 

(h(t)  decreases).  This  implies  that  the  unstable  manifolds  Wu(s^)  of  the 
saddles  of  period  k  intersect  the  stable  manifolds  Ws(sk  ^)  of  those  of 
period  k  -  1  and  one  therefore  obtains  the  global  structure  of  connected 
resonance  bands  indicated  in  Figure  C.3  (cf.  Hayashi  [1980],  Figure  3). 
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Figure  6.3  Heteroclinic  orbits  connecting  neighboring  resonance  bands 
m  =  3  and  m  =  2. 


This  structure  is  repeated  between  each  pair  of  resonance  bands  and,  using 

the  lambda  lemma,  it  is  possible  to  prove  that,  in  any  neighborhood  U  of 

the  stable  manifold  Ws(sk)  of  a  saddle  of  period  k  there  will  be  pieces 

of  the  stable  manifolds  Ws(s  . ),Ws(a . )  of  saddles  and  sinks  of  all  lower 

0  J 

periods  j  <  k.  The  resultant  violent  winding  and  packing  of  the  alternating 
domains  of  attraction  of  the  stable  subharmonics  leads  to  a  sensitive  depen¬ 
dence  on  initial  conditions,  especially  for  the  higher  order  subharmonic 
motions,  since  the  basic  width  of  the  domains  of  attraction,  calculated  from 
the  integral  ('6.4)  of  the  0(vD  Hamiltonian  approximation  is 

A! (m)  =  cos(m0)  -  iJ^m)^)^1^2  +  0(e)  (6.7) 


where  $  is  the  root  of 


60, (m) 

sin*  3  -  ^ys  -6R  (•*») 


(6.8) 


corresponding  to  the  center.  For  m  large,  while  J^(m),  JgCm.uj)  are  bounded 
(Theorem  3.8  and  §4),  fi'(m)  ■*  «  and  we  find 


Al Cm)  ~ 


(6.9) 


and  thus  the  domains  shrink  rapidly  in  size  as  m  increases.  In  typical 
computations  all  but  the  lowest  order  (m  <  5,  say)  subhamonics  lie  within 
a  band  whose  area  is  less  than  1%  of  the  interior  of  r+! 

In  Figure  6.4  we  illustrate  our  results  with  numerical  computations  of 
the  Poincare  map,  showing  unperturbed  resonant  orbits  and  the  perturbed  sub- 
harmonics.  For  more  detailed  and  extensive  computations  on  different  versions 
of  Duffing's  equation,  see  Hayashi  [1975,1980]  or  Ueda  [1980,  1981a].  The 
interaction  of  the  subharmonics  of  order  m  *  2  and  3  with  the  fixed 
points  near  (u,v)  =  (±1,0)  has  been  studied  using  conventional  averaging 
methods  by  Holmes  and  Holmes  [1981]. 

6.2  Homoclinic  Orbits  and  Chaotic  Motions 

As  we  pointed  out  in  Section  6.1,  the  results  based  on  the  averaged  sys¬ 
tems  (3.29)  and  (3.33)  are  invalid  near  the  homoclinic  orbit,  since 
|n'(m)|  -*■  «  as  m  <».  We  can,  however,  use  the  abstract  results  of  Section 
5  to  study  the  global  structure  of  the  Poincare  map  in  this  region,  and  we 
do  know  from  Theorems  3.6  and  3.8  and  Corollary  3.7  that  subharmonics  of  all 
orders  exist  here  for  y  >  R0(uj)<5,  although  we  cannot  compute  their  stability 
types  from  Eq.  (3.33). 

We  start  by  noting  that  the  symmetry  of  solutions  of  (4.1)  (if  (u,v,t) 
is  a  solution,  then  (-u,-v,t  +  tt/w)  is  also  a  solution)  together  with  the 
presence  of  transverse  homoclinic  points  for  y/6  >  R°(u>)  immediately  imply 
that  we  can  find  strips  S+,  S_  '  parallel*  to  Ws(p  )  and  in  a  neighborhood 


(b)  Subharmonics  of  orders  three  and  five  outside  H  *  0  and  of  orders  one 
and  three  inside  H  =  0. 

Figure  6.4.  Some  numerically  computed  subharmonics  for  Duffing' s  equation: 
to  3  4.2,  ey  =  0.65,  e6  3  0.005. 


Uv(p£)*  whose  images  P£(S+)  under  some  iterate  P£  of  the  Poincare  map 
are  as  shown  in  Figure  6.5.  The  orbit  of  any  point  x  e  S+  (resp.  S_)  will 
remain  near  the  unperturbed  loop  r°  (resp.  r°)  until  it  re-enters  IK 
In  this  way,  using  the  symbolic  dynamics  approach  of  Section  5.1,  we  can  find 
orbits  of  P  which  visit  neighborhoods  of  r°  and  r°  in  any  specified 
order. 


This  goes  some  way  toward  explaining  the  chaotic  snap  through  motions 
of  Figure  1.3  (cf.  Tseng  and  Dugundji  [1971],  Moon  and  Holmes  [1979],  Moon 
[1980]),  but  we  must  recall  that  all  these  orbits  are  of  unstable  saddle 
type,  their  dynamics  being  dominated  by  the  linearized  map  DP£(p£).  They 
would  thus  be  expected  to  give  rise  only  to  transient  chaotic  motions  and 
this  is,  indeed,  precisely  what  is  observed  when  numerical  integrations  are 
performed  in  certain  parameter  ranges;  almost  all  orbits  converging  to  fixed 
points  o»  subharmonics  lying  within  r°  or  (cf.  Figure  6.4).  However, 
in  view  of  Newhouse's  results  on  persistent  tangenices  and  wild  hyperbolic 
sets,  it  is  reasonable  to  expect  that  there  exist  stable  long  period  motions 
which  approximate  such  saddle-type  chaotic  orbits  arbitrarily  closely  on 
finite  time  intervals.  We  return  to  this  point  in  Section  6.3. 
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Note  that  while  such  'chaotic  subharmonics'  cannot  be  found  by  the  per¬ 
turbation  procedures  of  Section  3,  since  they  pass  through  Uv(p£ )  arbitrarily 
often,  their  existence  cajn  be  inferred  from  our  perturbation  calculations  of 
WS(P£)  and  WU(p£). 

We  also  note  that,  by  arguments  similar  to  those  of  Section  6.1,  we  can 
conclude  that  pieces  of  the  unstable  manifold  Wu(p£)  accumulate  on  all  the 
unstable  manifolds  Wu(sk),  k  =  1,  2,  —  of  all  periodic  saddle  type  mo¬ 
tions  within  r°  and  r°.  It  follows  that  any  attracting  periodic  motions 
within  this  region  lie  in  Ci,(Wu(p£)}. 

In  earlier  work  (Holmes  [1979])  we  showed  that  all  orbits  of  the  damped' > 

Duffing  equation  remain  bounded  and  enter  some  bounded  simply  connected  set 

1  2  i  t 

A  x  s  c  ir  x  s  as  t  ®.  Letting  A  e  Z  0  be  a  section  of  such  a  set, 

if  we  define  the  attracting  set  A  as 

A  =  n  (P  )n(A);  (6.10) 

n  >0  e 

it  follows  that 

A  =  C*(Wu(p  )}  .  (6.11) 

Here  assume  HfcCt  Orw  attyooh'il  A  twi"  oiusuv-  a.rf 

£ia  ctcSC  So  fpitm s  <r(-  Ue  U  fnfitM-pmxS* 

However,^  we  show  in  the  next  section,  it  is  not  likely  that  A  is  a  ‘nice1 

attractor;  an  invariant  attracting  set  which  is  indecomposible  and  contains  a 

dense  orbit.  Certainly  for  parameter  values  (y,6,u)  for  which  tangencies 

occur  Ct(Wu(p£)  contains  periodic  sinks  and  hence  is  decomposible. 


6.3  The  Strange  Attractor:  Numerical  Results  and  Conjectures 

We  start  by  summarizing  the  results  of  numerical  integrations  and  com¬ 
putations  of  the  Poincare  map,  drawing  in  the  results  of  Holmes  (1979],  Shaw 
[1980]  and,  principal ly,  Ueda  [1981b].  In  Figure  6.6a-d  we  show  plots  of 


the  stable  and  unstable  manifolds  Ws(p  ),  Wu(p  )  for  fixed  e6  *  0.25, 

C  w 

a)  *  1.0  and  ey  varying.  These  curves  are  obtained  by  iterating  short 
line  segments  lying  near  W*oc(p£),  W^0C(P£)*  We  note  that  the  manifolds 
have  just  intersected  for  the  first  time  at  ey  a  0.19  the  theoretical  value 
from  Eq.  (4.12)  being  ey  =  0.188.  As  y  is  increased  for  fixed  6  the 
homoclinic  bifurcation  and  resulting  formation  of  horseshoes  near  Wu(p£) 
occurs  while  stable  fixed  points  continue  to  exist  within  each  of  the  (unper¬ 
turbed)  loops  r®,  r®.  However,  with  appreciably  higher  y  and  for  a  wide 
range  of  y,6  values,  numerical  computations  suggest  that  there  are  no 

stable  low  period  attractors  and  one  observes  the  characteristic  'strange 
* 

attractor1  motions  of  Figure  6.7.  Comparing  Figure  6.7a  with  6.6d,  we  note 
that  this  single  orbit  appears  to  fall  on  CA  Wu(pe),  as  predicted  In  Section 
6.2.  For  typical  power  spectra  of  these  motions  see  Holmes  [1979]  or  Ueda 
[1980,  1981a]. 

In  Figure  6.8  we  indicate  regions  in  ey,  e<5  parameter  space  for  <d  ®  1 
in  which  the  various  attractors  are  observed.  We  also  superpose  some  bifur¬ 
cation  curves  from  our  theory.  In  each  region  we  sketch  a  representative 
motion  projected  onto  u-v  space.  These  numerical  results  are  due  to  Ueda 
[1981b],  Note  how  the  subharmonic  bifurcations  all  occur  very  close  to  the 
homoclinic  bifurcation  curve  y  =  R°(u>)6,  apart  from  those  of  order  m  »  1. 
Above  y  =  R^{m)6  subharmonics  of  all  orders  exist  both  inside  and  outside 
the  unperturbed  seperatrix  r*  u  {pQ}  u  f‘,  although  only  low  period  ones  are 

A  1 

observable.  The  theoretical  value  of  R  agrees  well  with  that  observed. 

■it 

This  term  is  used  only  for  convenience;  we  should  more  correctly  say  'appar¬ 
ent!.''  non-periodic  attracting  motion1.  As  indicated  at  the  end  of  Section 
6.2,  we  define  a  strange  attractor  for  the  map  P  as  an  attracting  set  A 
which  is  1 ndecomposible,  Has”  a  dense  orbit  and  is  neither  a  fixed  point  nor 
a  periodic  orbit. 


WW*-1"*1** '  1  m—"  —  I—.L MiiWi 


Ficure  6.7.  Numerically  computed  orbits  of  the  Poincare  map,  w  »  1.0, 
cy  =  0.3G.  (a)  =  0.25?  (b)  e6  -  0.20. 
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For  small  damping  (e6<  Q.1 ), sustained  strange  attractor  motions  are  not  observed 
and  we  find  only  fixed  points  (T-periodic  motions)  or  subharmonics  of  rela¬ 
tively  low  orders  (<  m  3  5)  coexisting  with  the  transient  chaos  characteris¬ 
tic  of  horseshoes,  cf.  Figure  6.4.  On  the  other  hand,  for  higher  damping 
levels  we  observe  stable  T-periodic  motions  outside  r*  u  {pQ}  u  r°  (cf.  Holmes 
[1979]).  One  apparently  requires  moderate  forcing  and  damping  for  the  strange 
attractor  to  exist.  We  note  that  in  some  parameter  ranges  (with  multiple 
shading  in  Fig.  6.8)  two  or  more  attractors  (periodic  and/or  strange)  coexist 
and  that  one  can  observe  hysteresis  in  the  jumps  from  one  attractor  to  another 
as  the  parameters  (e,y,6,oj)  are  slowly  varied. 

We  now  attempt  to  describe  the  various  attracting  sets  in  the  light  of 
the  theory  developed  in  this  article.  For  low  damping  the  stable  subharmonic 
,  motions  and  horseshoes  predicted  in  Sections  3  and  4  account  reasonably  well 

for  the  numerical  observations,  especially  in  view  of  the  small  characteristic 
•  size  of  the  domains  of  attraction  of  high  subharmonics  (eg.  (6.9))  and  the 

associated  implication  that  small  machine  errors  will  be  sufficient  to  per¬ 
turb  orbits  from  the  stable  manifold  of  one  such  subharmonic  to  a  neighboring 
one,  leading  to  a  'transitionally  chaotic  orbit*  (cf.  Ueda  (1980]).  Thus, 
although  there  certainly  are  parameter  values  for  which  such  long  period 
sinks  exist,  they  are  not  observable  in  practice.  Franks  [1981]  has  some 
interesting  results  on  the  addition  of  random  perturbations  to  a  motion 
originally  attracted  to  a  stable  oeriod  7  sink  of  the  Hdnon  map  (cf.  Hdnon 
[1976]).  The  addition  of  a  small  random  vector  after  each  iteration  leads 
to  an  orbit  which  appears  to  fill  out  the  same  ’strange  attractor'  as  one 

( 

obser/<?s  for  nearby  parameter  values  without  intentional  perturbations.  This 
observation  is  one  of  the  motivations  for  tne  conjecture  which  follows. 
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The  second  pertinent  observation  is  also  based  jointly  on  analysis  and 
numerical  integrations.  The  study  of  near  resonant  excitation  of  the  stable 
periodic  motions  near  (u,v)  =  (±1 ,0)  in  Holmes  [1979]  and  more  particularly 
Holmes  and  Holmes  [1981]  shows  that  for  forcing  frequencies  <u  near  or 
2/7,  one  obtains  harmonic  and  subharnrnic  bifurcations  as  y  is  increased, 
so  that  the  stable  fixed  point  of  the  Poincard  map  can  become  a  saddle  point 
and  (in  the  second  case)  throw  off  an  orbit  of  period  two.  As  in  the  case 
of  subharmonics  near  the  homoclinic  orbit  (cf.  §5.2),  there  is  some  reason 
to  think  that  this  period  doubling  bifurcation  is  the  first  in  a  countable 
sequence  which  has  an  accumulation  point  (cf.  Feigenbaum  [1978]  Collet, 
Eckmann  and  Lanford  [1980]).  The  studies  of  Equation  (4.1)  and  an  associated 
cubic  mapping  of  the  plane 

(u,v)  -*•  (v,  -  8u  +  av  -  v3)  (6.12) 

strongly  support  this  idea  (  Holmes  [1979]).  In  fact  in  view  of  the  results  of 
Gavrilov  and  Silnikov  outlined  in  Section  5.2,  we  can  conclude  that  aVl_  the 
subharmonic  sinks  of  sufficiently  high  period  2um/tu  undergo  such  period 
doubling  bifurcations  as  y  increases  for  fixed  6.  It  seems  reasonable  to 
expect  these  each  tc  be  the  first  of  such  an  accumulating  sequence,  and  all 
such  sequences  presumably  accumulate  on  the  homoclinic  bifurcation.  This 
provides  a  mechanism  by  which  all  sinks  of  periods  less  than  some  fixed 
integer  N  can  become  unstable,  leading  to  our  conjecture: 

Conjecture  6.1.  For  any  fixed  Integer  N  <  there  are  open  sets  of  param¬ 
eter  values  in  (c.y.Oiw)  space  for  which  the  attracting  set  A  a  Pn(A) 
Cl(Hu(p  ))  of  the  Duffing  equation  (4.1)  contains  a  shift  on  two  symbols 
with  a  countable  set  of  saddle  type  periodic  orbits  and  finite  sets  of 


periodic  sinks,  none  of  the  latter's  periods  being  less  than  N. 
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This  conjecture  implies  that  A  is  not  a  strange  attractor,  since  It  can  be 
decomposed  into  a  finite  set  of  simple  periodic  attractors.  Such  an  attract¬ 
ing  set  is  structurally  stable.  We  can  of  course,  prove  this  conjecture  for 
the  case  N  =  2  (Holmes  and  Holmes  [1981]).  The  conjecture  can  be  proven 
for  arbitrarily  high  N  using  a  modified  Duffing-van  der  Pol  equation  (Holmes 
[1981]).  Numerical  work  of  Chirikov  [1979,1980]  and  Lieberman  [1981]  on 
diffeomorphisms  of  the  plane,  also  lends  support  to  our  conjecture. 

Newhouse's  results  [1974,  1979,  1980]  suggest  that  there  are  also 
residual  subsets  of  parameter  values  for  which  one  has  countably  many  sinks, 
but  a  careful  examination  of  the  proof  of  his  Theorems  5.4  and  5.5  shows 
that  the  small  perturbations  of  diffeomorphisms  with  quadratic  tangencies 
necessary  to  obtain  such  infinite  sets  of  sinks  are  very  special,  and  there 
appears  to  be  no  guarantee  that  our  specific  family  realizes  these  perturba¬ 
tions.  In  fact  some  recent  work  of  Guckenheimer  [1981]  which  borrows  ideas 
for  Yakobson's  [1978  ]  work  on  one  dimensional  maps  indicates  that  one  can 
also  construct  diffeomorphisms  close  to  such  homoclinic  bifurcations  which 
possess  no  periodic  sinks  and  therefore  presumably  have  a  genuine  (indecom¬ 
posable)  strange  attractor.  (Such  an  attractor  has  been  proven  by  Misiurewicz 
[1980]  to  exist  for  the  piecewise  linearLozi.  map.  Here  the  persistent  quad¬ 
ratic  tangencies  and  their  associated  sinks  do  not  occur  since  the  folds  in 
the  manifolds  become  points  in  the  piecewise  linear  manifolds).  On  the  basis 
of  this  work  one  can  conjecture  that  a  situation  similar  to  that  in  one  dimen¬ 
sional  quadratic  maps  holds: 

Conjecture  6.2.  There  is  a  set  of  parameter  values  of  positive  measure  but 
containing  no  open  sets  for  which  the  Duffing  equation  has  a  strange  attractor 
A = Cl (Wu(p£) ) ,  i.e.  a  contains  a  dense  orbit  and  is  indecomposable. 

Note  that  our  two  conjectures  are  compatible. 

On  this  note  of  unbridled  speculation,  let  us  close  our  story. 


APPENDIX.  INVARIANT  MANIFOLDS  AND  THE  LAMBDA  LEMMA 
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Stable  Manifold  Theorem.  Hirsh  Pugh -Shub  [1977].  Suppose  that  p  is  a 
hyperbolic  fixed  point  of  a  Cr  di ffeomarphism  P  :  IRn  -*■  Kn  and  that 
Es (p) ,  Eu(p)  are  the  stable  and  unstable  aigenspaces  of  the  associated 
linearized  mapping  DP ( p ) .  Then,  in  a  neighborhood  U.(p)  there  are  Cr 
invariant  submanifolds  '^oc(p)>  wifoc(p),  tangent  to  Es(p),  Eu(p)  at.  p. 
The  global  manifolds  Ws(p),  Wu(p)  are  injectively  immersed  copies  of 
IRS ,  IRU  ,  when.  s  =  dim  ES(p),  u=dim  Eu(p)  (s  +  u  =  n).  Moreover,  if  P 
depends  smoothly  upon  parameters  p  £  IRm,  so  do  w*oc(p)>  W^tp). 

X- Lemma .  (Palis  [1969]).  Let  P  be  a  Cr  diffeomorphism  with  a  hyperbolic 
fixed  point  p  and  let  Du  be  a  u-disc  in  Wu(p).  Let  A  be  a  u-disc 
meeting  WS(p)  transversely  at  some  point  q.  Then  U  Pn(A)  contains 
U-discs  arbitrarily  Cr  close  to  Du. 

These  results  also  apply  to  periodic  orbits  (p.  }i  =1,  ....  n 
Pn(Pi)  a  Pi  and  an  analogue  of  the  first  result  applies  to  flows,  such  as 
that  of  the  suspended  system  (2.2). 
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